Implications of a 300 — 500 GeV/c 2 Z' boson on pp collider data at ^fs = 1.8 TeV 
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Recent analyses of precision low energy electroweak data indicate that the deviation from the 
Standard Model predictions of the measurement of atomic parity violation (2.3 a), the effective 
number of massless neutrinos (2 a), and Ab (2.7 a) could be better described if the existence of an 
extra Z' neutral gauge boson is assumed. We investigate the implications of a 300 - 500 GeV/c 2 extra 
Z' on current pp collider data at \/s = 1.8 TeV, including the forward-backward charge asymmetry 
for very high mass e + e~ pairs, and the invariant mass spectrum of high mass e + e~ , p + \i~ , ti and bb 
pairs. For example, a 500 GeV/c 2 Z' with a total production cross section of « 3pb and enhanced 
coupling to the third generation, better describes both the low energy and the Tevatron data. 
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Recent analyses of precision electroweak data in- 
dicate that there are several measurements for which 
the deviations from the Standard Model predictions are 
larger than two standard deviations (a). These include 
the measurements of atomic parity violation a\ (2.3 a), 
the effective number of massless neutrinos J| (2 a), and 
Af, j|] (2.7 a). These analyses show that the data are 
better described if an extra Z' neutral gauge boson is 
assumed. In this analysis, we investigate if there is 
evidence for a Z 1 boson in current pp collider data at 
y/s = 1.8 TeV. 

We take the parameters and couplings of the Z' to 
the first generation quarks and leptons from the analysis 
of the atomic parity violation data by Rosner Q . Ros- 
ner's analysis indicates that the atomic parity violation 
data H are better described with an Eg extra Z' bo- 
son. Within this model || there is a continuum of Z' 
possibilities given by Z' = Z^ cos0 + Z x smqb. Rosner 's 
analysis of the atomic parity violation data yields a best 
fit for a region of allowed Z 1 mass, Mz<, and cf>. For ex- 
ample, For cf> — 120° the data are best fit with an Eg Z' 
with a mass of about 800 GeV/c 2 . For an E 6 Z' with 
Mz> = 500 GeV/ c 2 , the atomic parity violation data are 
best fit with cf) = 70° ± 5° and = 160° ± 5°, and for 
a Z' with a mass of 350 GeV/c 2 , these data are best fit 
with (j) = 60° ± 5° and (j) = 173° ± 5°. 

Erler and Langacker H extend the Z' analysis to in- 
clude all precision electroweak data, and include more 
general classes of Z' models. In one of the cases, the 
analysis is extended to allow for the coupling to the third 
generation to be different from the coupling to the first 
two generations. With these additional parameters, they 
are not able to place a constraint on Mz>, but the low 
energy electroweak data prefer a Z' with a small (but 
finite) mixing to the Z , and a large coupling to the third 
generation, as expected in some models ^||. The larger 
coupling to bb pairs is needed to account for the 2.7 a 
deviation of Q| from the Standard Model prediction. 

Although the mass limits @ from CDF and DO || for 
a variety of Z' models are in the 600 GeV/c 2 range, the 



limits are reduced § by 100 to 150 GeV/c 2 , if the Z' 
width (typically Yz> ~ 0.01 Mz') is increased to account 
for the possibility of additional decays modes to exotic 
fcrmions (which are predicted in Eq models), and/or su- 
persymmetric particles. The limits are even lower if one 
includes the possibility of a more general model with en- 
hanced couplings to the third generation. Therefore, we 
investigate the present Run I collider data for high mass 
e + e~, it, and bb final states, and look for pos- 

sible signatures for a Z' extra gauge boson of the kind 
that is favored by the low energy data. We constrain the 
relationship between the couplings to the first two gen- 
erations and the mass of the Z' to be the same as that 
for an Eq Z' boson from Rosner's fits to the low energy 
measurements. 

In hadron-hadron collisions at high energies, massive 
e + e~ and pairs are produced via the Drell-Yan 

7* jZ process. The presence of both vector and axial- 
vector couplings in this process gives rise to an asymme- 
try Q, Afb, in the final-state angle of the lepton in the 
rest frame of the e + e~ and /i + pT pair (with respect to 
the proton direction). Within the Standard Model, for 
M ^> Mz, the large predicted asymmetry (w 0.61) is 
a consequence of the interference between the propaga- 
tors of the 7* and Z. New interactions such as an extra 
Z' boson result in deviations from the standard model 
predictions in both da/dM and Afb- 

Figures [l] and |^ compare the measured high mass 
Drell-Yan da/dM (CDF §|l0) and D0 [§) and A FB 
(CDF |lCf| ) to theoretical predictions. The Standard 
Model da/dM curve is a QCD NNLO calculation 
with MRST99 NLO PDFs @. The predictions in Fig- 
ure |l](a) are normalized by a factor F = 1.11, the ra- 
tio of the CDF measured total cross section in the Z 
region |l3| to the NNLO prediction (the overall normal- 
ization uncertainties arc 3.9% for the experimental data 
and 5% for the NNLO theory). The Standard Model pre- 
diction for Afb has been calculated (TJ] in QCD-NLO 
with sin 2 #gjPj = 0.232. The measured da/dM and Afb 
values are in good agreement with the standard model 



1 



predictions. In the two highest mass bins (4 events in 
the 300 - 600 GeV/c 2 range), A FB is about 2.2 standard 
deviations below the standard model prediction (there 
are 3 events in the negative hemipshere and one event 
in the positive hemisphere). An asymmetry in the 300 - 
600 GeV/c 2 range which is smaller than the Standard 
Model prediction could result from the exchange of a 
300 - 500 GeV/c 2 Z' gauge boson. 
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FIG. 1. (a) da/dM distribution of e+e" (CDF and D0) 
and fi + n~ pairs (CDF). The Standard Model theoretical pre- 
dictions (dashed line) have been normalized (by a factor of 
1.11) to the CDF data in the Z boson mass region, (b) CDF 
Afb versus mass compared to the standard model expecta- 
tion (dashed) . Also shown are the predicted theoretical curves 
(xl.ll) for da/dM and Afb with an extra E@ boson with 
M z > = 350 GeV/c 2 and V z , = 0.1 M z > , for 4> = 60° (solid) 
and 4> = 173° (dotted). 

For the Eq Z' models that were fit to the low energy 
electroweak data, the couplings to the first two gener- 
ations are well constrained. Contributions to the Z' 
cross section from tt and bb annihilation are strongly sup- 
pressed due to the small it and bb parton luminosities, 
even in the case where the Z'ti and Z'bb couplings are 
strongly enhanced. For a given Z' mass and <f>, the total 
Z' production cross section (in all channels) in pp col- 
lider data at y/s = 1.8 TeV thus is determined by the 
couplings of the Eq Z 1 to the first and second generation 
up and down type quarks. The partial width to elec- 
trons is also determined. The integrated cross section 
for e + e~ final states, <txBR(Z' — * e + e~), is determined 
by the Z' — > e + e~ branching ratio, and is therefore pro- 
portional to 1/Tz'- In contrast, the prediction for Apg, 
which results from the interference between the Standard 



Model and the Z' amplitudes, is quite insensitive to the 
Z' width. In this study, we compare the Run I collider 
data to a model with a Z' width of Tz< =0.1 Mz> (which 
is about a factor of 10 larger than the expected width for 
an Eq Z' boson for the case of universal couplings to all 
three generations) . This allows for enhanced couplings to 
the quarks of the third generation, or for additional decay 
modes to exotic fermions or supersymmetric particles. 
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FIG. 2. Same as Figure but shown here are the predicted 
theoretical curves for da/dM (xl.ll) and Afb with an extra 
E 6 boson with M z , = 500 GeV/c 2 and Y z , = 0.1 M Z ', for 
(f> = 70° (solid) and 4> = 160° (dotted). 

The predictions for da/dM and Afb for the case of the 
Drell-Yan process including an extra Eg Z' boson (for a 
given M and <j)) are first calculated in leading order in 
QCD (with MRSR2 NLO PDFs), ignoring contributions 
from tt and bb annihilation. da/dM is then multiplied 
by a differential if -factor, K(M), which is determined 
by comparing the leading order prediction for da/dM 
including a 50 TeV/c 2 Z' (i.e. effectively no Z' in the 
model) to the Standard Model da/dM prediction from 
the QCD NNLO 0] calculation (with MRST99 NLO 
PDFs). In the 300 - 400 and 400 - 600 GeV/c 2 ranges, 
we obtain a differential if-factor of 1.172 and 1.154, re- 
spectively. Since we plan to compare to CDF Drell-Yan 
data, the predictions in Figure |(a) and|(a) are normal- 
ized by a factor F = 1.11, the ratio of the CDF measured 
total cross section in the Z region [jl3| to the NNLO pre- 
diction. Because the data show a clustering of events at 
mass values of 350 and 500 GeV/c 2 , we investigate the 
signatures for a Z' boson with these two mass values. As 
mentioned earlier, in all of the calculations we use a Z' 
width of Y z> = OAMz'- If we use a smaller width, the 
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cross section in the di-lepton channel would be larger by 
a factor proportional to 1 /Y Z ' ■ 

For an E & gauge boson with M Z i — 350 GeV/c 2 , 
IV = O.lMy, and <p = 60° (173°), the theoretical pre- 
diction (xl.ll) for the integrated total cross section of 
e + e~ pairs in the 300 — 400 GeV/c 2 mass range is 73 
(31) fb. The corresponding cross section for the Stan- 
dard Model prediction (xl.ll) in this range is 54 fb. 
Both the Standard Model cross section, and the cross 
section including an additional 350 GeV/c 2 Z 1 are con- 
sistent with the observed CDF cross section in this range 
of 81±47 fb. Our predicted theoretical curves for da/dM 

and Afb with an extra Eq boson, for Mz* = 350 GeV/c 2 a Z' boson with 500 GeV/c 2 varies from 3.1 pb for a very 



r 'z' = 0.1 Mz>, are shown in Figure | The total (with 
decay to all channels) production cross section (xl.ll) 
for this Z 1 , together with the Z' — > e + e~ branching ratio 
and the Standard Model Drell-Yan cross section in the 
400 — 600 GeV/c 2 mass range are listed in Table || For 
Afb, the prediction for = 70° gives better agreement 
with the data than that for <f) = 160°. For <j) = 70° (160°), 
the probability that theory and data agree is 6.3% (3.7%) 
(see Table |). 

We now look for possible signatures of a Z' with these 
production cross sections in the bb and ti channels. In 
the bb channel, the CDF 95% CL cross section limit for 



and Tz> = O.lMy, are shown in Figure H The total 
production cross section (xl.ll) for this Z' for <fi = 60° 
(173°) is 17 (11) pb. For A FB , <t> = 60° results in a bet- 
ter agreement of theory and data than <fi = 173°. In 
the 300 — 600 GeV/c 2 range, the probabilty that the for- 
ward backward asymmetry predicted by the Standard 
Model agrees with CDF data is 2.4% (see Table 1. For 
M Z ' = 350 GeV/c 2 , T z > = 0.1 M z >, and <j> = 60° (173°), 
the corresponding probability is 13.4% (7.7%). 

Because the couplings of the E§ Z 1 to the first two 
generations of quarks and leptons are constrained by the 
fit to the low energy electroweak data, most of this cross 
section should appear in the form of decay modes to ei- 
ther exotic fermions (which are predicted in Eq mod- 
els) and/or supersymmetric particles, or decays to third 
generation quarks. Since the analysis of Erler and Lan- 
gacker [|J indicates that a larger coupling to the third 
generation is needed to account for the 2.7 a deviation 
from the Standard Model of At, j| , we look for Z' signa- 
tures in the invariant mass spectra of ti and bb high mass 
pairs. The CDF 95% CL limit on the bb cross section 
for a 350 GeV/c 2 Z' varies from 12 pb for a very narrow 
Z' to 28 pb for a Z' with T z < = 0.3 M z >- Therefore, a 
production cross section of 11 to 17 pb for a 350 GeV/c 2 
Z' with Tz' = 0.1 Alz' (which predominantly decays to 
bb pairs) is consistent with these bb limits. Note that 
Z' — > ti decays are either forbidden or strongly phase 
space suppressed for Mz> — 350 GeV/c 2 . 



narrow Z' to 5.5 pb for a Z' with Tz> — 0.3 Mz>- In 
the ti channel, the CDF 95% CL cross section limit for 
a Z' boson with M z > = 500 GeV/c 2 is 7.5 pb. There- 
fore, a production cross section of 1.7 to 3.2 pb for a 
500 GeV/c 2 Z' with T z . =_0AM z < (see Table |0j), which 
predominantly decays to bb and/or ti pairs, is consistent 
with these limits. It is interesting to note that the pub- 
lished ti and bb mass distributions show a slight excess 
of events in the 500 GeV/ c 2 region. 
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FIG. 3. The invariant mass distribution of tt pairs at 
CDF (fit with m t = 175 GeV/c 2 ). The published CDF Z' 
cross section 95% CL upper limit in the tt channel is 7.5 pb 



■ 2 mass range and the Z' — > e + e 
= 350 GeV/c 2 and T z > = 0.1 M z 



The Standard Model e 
the 300 - 400 GeV/c 
cross section for M z < 

with (j) = 60° and (f> = 173° are summarized in Table |J. 
The Z' — > e + e~ branching ratio and the total Z' produc- 
tion cross section are also shown in the table. 

For an E 6 gauge boson with M z < = 500 GeV/c 2 , 
Y z > = 0.1 M z >, and <j> = 70° (160°), the theoretical pre- 
diction (xl.ll) for the integrated total cross section of 
e + e~ pairs in the 400 — 600 GeV/c 2 mass range is 17 
(7) fb. The corresponding cross section for the Standard 
Model prediction (xl.ll) in this range is 17 fb. Both the 
Standard Model cross section, and the cross section in- 
cluding an additional 500 GeV jc 2 Z' are consistent with 
the observed CDF cross section in this range of 28±28 fb. 
Our predicted theoretical curves for da/dM and Afb 
with an extra Eq boson, for M z > = 500 GeV/c 2 and 



Drell-Yan cross section for at M it = 500 GeV/c . The 1.4 a excess at 500 GeV/c 2 in 

the invariant mass spectrum corresponds to a cross section of 
2.3±1.7p6. 



Figure || shows the CDF published [jl6],[l7| invari- 
ant mass distribution of ti pairs (assuming m t = 
175 GeV/c 2 ). The 1.4 a excess at 500 GeV/c 2 in the in- 
variant mass spectrum corresponds |l7]] to crxBR(tt) = 
2.3 ± 1.7 pb. Therefore, the CDF ti data support the 
hypothesis of a Z' boson with a large coupling to the 
third generation. Although DO has not searched for res- 
onances in the ti channel, the published D0 mass spec- 
trum JlJ] for ti events is consistent with an enhance- 
ment in the 460 — 500 GeV jc 2 region at a similar level. 
There is also small 1 a excess at 500 GeV jc 2 in the 
CDF bb invariant mass spectrum which corresponds to 
crxBR(66) = 1 ± 1 pb. Therefore, the CDF bb data are 



3 



also consistent with the hypothesis of a Z' boson with a 
total production cross section of 1.7 to 3.2 pb which has 
a small branching ratio to di-leptons and predominantly 
decays to top and bottom quarks. Note, that for a Z' 
that mixes with the Z, the fits to low energy electroweak 
data at the Z peak already constrain || the level of the 
coupling of a Z' to r leptons to be similar to the cou- 
plings to electrons and muons. Therefore, an enhanced 
signal in the t + t" channel is not expected for a Z' which 
mixes with the Z. 

In summary, we find that cither a 350 or a 500 GeV/ c 2 
extra Z' with Yz> = 0.1 Mz' and enhanced couplings 
to the quarks in the third generation not only gives a 
better description of the low energy electroweak data, 
but also better describes the forward-backward asymme- 
try for e + e~ pairs in the 300-600 GeV/c 2 range. A 
500 GeV/c 2 extra Z' also accounts for the 1.4 a excess 
at 500 GeV/c 2 in the invariant mass spectrum of high 
mass of tt events at CDF. A Z' in the 350 - 500 GeV/c 



mass range is also compatible 19 with the latest mea- 
surement of the muon anomalous magnetic moment p0[ . 
With the upgraded CDF and D0 detectors, and the antic- 
ipated factor of 20 higher luminosity in Run II, improved 
searches for Z' bosons could be made in the invariant 
mass spectrum and in the forward-backward asymme- 
try for all three di-lepton channels (ee, mi and tt). In 
addition, the improved silicon vertex detectors in both 
experiments will increase the sensitivity of such searches 
in tt and bb final states. 
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TABLE I. The measured CDF (preliminary) and predicted 
average Afb in the 300 — 600 GeV/c 2 range and the proba- 
bility that models can result in the observed Afb (3 events 
in the negative hemisphere and 1 event in the positive hemi- 
sphere) for the Standard Model and models with an extra Z' 
boson. 



M z ,{GeV/c 2 ) 




A FB 


Probability 


data 




-0.45 ±0.47 




SM 




0.612 


2.4% 


350 


60° 


0.284 


13.4% 


350 


173° 


0.419 


7.7% 


500 


70° 


0.458 


6.3% 


500 


160° 


0.553 


3.7% 
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TABLE II. CDF (Preliminary) data versus Standard 
Model Drell-Yan cross section a ee (fb) for two mass bins: 
300 - 400 and 400 - 600 GeV/c 2 . Also shown are total 
Drell-Yan cross sections when an additional Z' boson with 
a mass M z , with a total width of F z i — 0.1 M z i is included. 
The Z' di-lepton branching ratio BR(ee) and total production 
cross section a to t (pb) (with final state decays to any parti- 
cles) are calculated for the values of <j> obtained from Rosner's 
fits to low energy data. 



Bin (GeV/c 2 ) M z ,{GeV/c 2 ) 


4> 


a ee (fb) BR(ee) 


o-tot(pb) 


300 - 400 


data 




81 ± 47 




300 - 400 


SM 




54 




300 - 400 


350 


60° 


127 0.43% 


17.1 


300 - 400 


350 


173° 


85 0.29% 


10.9 


400 - 600 


data 




28 ± 28 




400 - 600 


SM 




17 




400 - 600 


500 


70° 


34 0.53% 


3.2 


400 - 600 


500 


160° 


24 0.41% 


1.7 
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